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ABSTRACT 

We present models of giant planet migration in evolving protoplanetary disks. Our disks evolve subject to 
viscous transport of angular momentum and photoevaporation, while planets undergo Type II migration. We 
use a Monte Carlo approach, running large numbers of models with a range in initial conditions. We find 
that relatively simple models can reproduce both the observed radial distribution of extra-solar giant planets, 
and the lifetimes and accretion histories of protoplanetary disks. The use of state-of-the-art photoevaporation 
models results in a degree of coupling between planet formation and disk clearing, which has not been found 
previously. Some accretion across planetary orbits is necessary if planets are to survive at radii < 1 .5AU, and if 
planets of Jupiter mass or greater are to survive in our models they must be able to form at late times, when the 
disk surface density in the formation region is low. Our model forms two different types of "transitional" disks, 
embedded planets and clearing disks, which show markedly different properties. We find that the observable 
properties of these systems are broadly consistent with current observations, and highlight useful observational 
diagnostics. We predict that young transition disks are more likely to contain embedded giant planets, while 
older transition disks are more likely to be undergoing disk clearing. 

Subject headings: accretion, accretion disks - planetary systems: formation - planetary systems: protoplane- 
tary disks 



1. INTRODUCTION 

Understanding how planets form has been an active topic of 
research for centuries, but interest in this subject has increased 
dra matically since the discovery of the first extra- solar plan- 
ets dMavor & Ouelozlll995l;lMarcv & ButleJl 19961) . Over 350 
such plan ets are now know n, with a diverse range of proper- 
ties (e.g., |UdrxitiDllQQ3). It was recognized very quickly 
that many extra-solar planets orbit very close to their parent 
stars, and that such planets could not have formed at their cur- 
rent locations. This in turn sparked re newed interest in the es- 
tablished theory of plane t migration (iGoldreich & Tremaind 
119801; iLin & Papaloizoul [1986). in which planets form far 
from their parent stars and "migrate" to smaller radii through 
the action of tidal torques. The migration of low-mass planets 
remains controversial, but so-called Type II migration, which 
applies to planets sufficiently massive (> O.SMjup) to open 
gaps in their parent gas disks, is now relatively wel l under- 
stood (see, e.g., the review by Papaloizou et al. 2007). 

In a similar vein, it was discovered in the 1980s that young 
stars, such as the T Tauri stars (TTs) , possess circumstellar 
disks (e.g. JSargent & Beckwithl[l987h . These disks are of or- 
der a percent of the mass of their central star dBeckwith et alj 
[rggO*; And rews & Williamsll2005h . and we now have a well- 
established evolutionary picture where disk-bearing, classi- 
cal T Tauri stars (CTTs) evolve into disk-less, weak-lined 
T Tauri stars (WTTs ) on timescales of a few Myr (e.g., 
iHartmann etal] 119981: iHaisch et alj l2001h . The dominant 
processes driving (gas) disk evolution are angular momen- 
tum transport , and evaporation due to heating by energetic 
photons (e.g.. lHollenbach et allllOOOl; iDullemond et alj|2007l: 
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lAlexande3l2008h . and these evolving "protoplanetary" disks 
are the sites of planet formation. 

Statistically, it is clear that some fraction of observed proto- 
planetary disks must contain planets. Identifying these planet- 
bearing disks is of considerable interest, since doing so would 
provide evidence as to where and when planets typically form 
within disks, and discriminate between different models for 
giant planet formation. Unfortunately, although we have a 
reasonable idea of what a disk containing a planet would look 
like (an ordinary disk at large radii, but depleted of gas and 
dust interior to the planet's orbit), we cannot say that all the 
disks observe d to show such signatures (the so-cal led tran- 
sitional disks: IStrom et alJll989E iNaiita et aT]|2007h contain 
planets. A significant fraction of transitional disks may in- 
stead represent an interme diate stage of disk evolution prior 
to final disk clearing (e.g.. lCieza et a l. 2008), a process that 
occurs in the presence or absence of planets. Here, we seek to 
construct models that include both disk evolution and planet 
formation. By comparing the models against both exoplanet 
and disk evolution statistics we seek to make maximum use 
of available observational constraints, and thereby predict the 
observational appearance of planets within evolving disks. 

In this paper we present models of giant planet migration in 
evolving protoplanetary disks. We restrict ourselves to con- 
sidering relatively massive planets, > O.SMjup, primarily be- 
cause searches for exoplanets are pr esently only complete to 
around the Jupiter mass level (e.g.. iFischer & Valentil 120051 : 
lUdrv et akl 12007) . Consequently, throughout this paper the 
term "planet" always refers to gas giant planets; we make 
no attempt to model planets of lower mass. Our models in- 
clude viscous transport of angular momentum, photoevapora- 
tion, and Type II planet migration. We adopt a Monte Carlo 
approach, running large numbers of models with a range of 
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initial conditions in order to follow the time evolution of dis- 
tributions of disk and planet properties. In ^we present our 
numerical model, and compare the results to the observed ra- 
dial distribution of extra-solar planets, and to a large range 
of observations of protoplanetary disks. Our results compare 
favourably with the observed properties of both planets and 
disks, and we discuss how the observational data can con- 
strain various properties of the model. This is the first study 
to model populations of transitional disks theoretically, and 
in ^ we discuss the transition disk phenomenon. Our rela- 
tively simple model reproduces the known properties of these 
objects well. The model produces transition disks via two 
different mechanisms (gap-opening by planets and disk clear- 
ing), and we discuss the relative efficiency of these processes 
in the models. We show that observations of transition disk 
masses and accretion rates remain the most straightforward 
means of distinguishing between different types of transition 
disks, and make predictions for future observations of such 
objects. 

2. MODELS 

2.1. Planet migration model 

In our model, protoplanetary disks evolve due to viscous 
transport of angular momentum and photoevaporation by the 
central star Planets migrate due to tidal interaction with the 
disk (in the Type 11 migration regime), and the disk is also 
subject to tidal torques from planets. The coupled evolution of 
a prot oplanetary disk and a pla net is described by the equation 
(e.g.. iLin & Papaloizoull 19861) 



dt 
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(1) 

Here S(/?,0 is the disk surface density, t is time, R is cylindri- 
cal radius, i/ is the kinematic viscosity, and = IM© is the 
stellar mass. The first term on the right -hand side describes or- 
dinary viscous evolution of the disk (iLvnden-Bell & Pringld 
ri974t IPringld 119811) . and the Ew(/?,f) term represents the 
mass-loss due to photoevaporation. The second term de- 
scribes how the disk responds to the planetary torque: here 
A{R, a) is the rate of specific angular momentum transfer from 
the planet to the disk . Following Trilling et al. (1998) and 
lArmita^e et al.l (12002 ). for a planet of mass Mp = qM^ at ra- 
dius (semi-major axis) a we adopt the following form for A 
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where 



Ap = max(H,\R-a\) 
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and H is the disk scale-height. This form for A is the same 
as that used by Lin & Papaloizou ( 1986), but modified to give 
a symmetric treatment inside and outside the planet's orbit. 
This transfer of angular momentum causes the planet to mi- 
grate at a rate 
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This treatment of planet migration is necessarily idealized, but 
has previously been shown to give res ults comparable to mo re 
sophisticated numerical models (e.g.. lTakeuchi et aLlll996h . 



The kinematic viscosity governs the transport of angular 
momentum in the disk, and we adopt an alpha-disk model 

i^(R) = aflH^ , (5) 

where a is the standard [Shakura & SunvaevI (Il973h viscosity 
parameter and ft(R) = yj GM^jR} is the orbital frequency. We 
adopt a scale-height consis tent with a flared disk model (e.g., 
iKenyon & Hartmannlll987h 

H(R) (X RP , (6) 

where the power-law index /? = 5/4. This choice of p gives 
a viscosity i/ oc R, consistent with high-resolution obser- 
vations of disk structure (e.g., Andrews & Williams '2007t 
Andrews et al. 2009). As in Alexander & Armitage (2007), 
we normalize this relationship so that the disk aspect ratio 
H/R = 0.0333 atR = lAU. This parametrization assumes that 
stellar irradiation dominates over viscous heating, and is thus 
strictly valid only for accretion rates less than a few times 
lO"**M0yr"' (e.g., D'Alessio et al. 1999). In practice, this 
means that our model of angular momentum transport in the 
disk is only physical at times t > lO^yr 
The rate of mass-loss due to photoevaporation is 
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The "diffuse" profile applies when the inner disk is opti- 
cal ly thick to ionizi ng photons, and was studied in detail 
by Hollenb ach et alj j 19941) and IFont et alj (12004 , see also 
iLiffman 2001 7" In this case radiation from the star creates 
a thin ionized layer on the disk surface, with a sound speed 
Cj ~ 10km s"', and the diffuse (recombination) radiation field 
from the disk atmosphere is the dominant source of ionizing 
photons at ~AU radii. Outside some critical radius the heated 
layer is unbound and flows as a wind. The mass-loss profile 
is strongly concentrated close to the critical radius 
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where = GM^^/cl is the "gravitational radius" defined by 
iHoUenbach et al.l (Il994t) ^. The "direct" profile applies when 
the inner disk has been cleared (by either viscous or tidal 
torques) and is optically thin to ionizing photons. In this case 
stellar irradiation ionizes the inner edge of t he disk directly, 
and the disk is cleared from the inside out ( Alexander et ajj 
|^06a.b). We make use of the numerical parametriza- 
tions given in the appendix of I Alexander & Armitagd (l2007l) . 
and switch between the two profiles when the surface den- 
sity in the inner disk {R < Ra-it) falls below some criti- 
cal value Sc. We set Sc = 10~^g cm"^, but note that the 
results are not ver y sensitive to the exact val ue adopted 
(see discussion in [Alexander & Armitagd l2007h . We as- 
sume a stellar ionizing fl ux of < 1> = lO'^^photons per sec- 
ond (e.g., [A lexander et al] l2005t iPascucci & SterzikI 120091: 
iHollenbach & Gortiii2009i) . which results in a (diffuse) wind 
rate of Mw 4 x lO-i°M0yr-i. 

The initial surface density profile is taken from the similar- 
ity solution of the diffusion equation dLvnden-Bell & Pringld 
IT974tllfotmann et alJll998h : 



E(R): 
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exp(-R/R,), 



(9) 



^ Sg is found by equating tlie Keplerian orbital speed with the sound speed 
of the ionized gas, Cs. 
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where Md is the initial disk mass. The scaling radius 
defines the initial disk size, and also sets the viscou s time- 
scale of the disk. Following i Hartmann et al.l (Il998h we set 
Rs = lOAU, which results in a viscous time-scale of ti, = 
Rl/3v(Rs) ~ 5 X lOV for a = 0.01. 



2.2. Numerical method 

We solve the diffusion equation for the gas surface den- 
sity (Equation [U numerically using a standard first-order ex- 
plicit finite difference scheme, on a grid of points equispaced 
in R^l^ (e.g.. lPringle et al.||1986h . We choose a grid covering 
the range [0.04AU, lOOOOAU], and adopt zero-torque bound- 
ary conditions (i.e., we set E = in the boundary cells). In 
the absence of a planet the gas disk can be evolved on a fairly 
coarse grid, but higher resolution is required to maintain ac- 
curacy when calculating the effect of the planetary torque. In 
order to make efficient use of computing resources we there- 
fore use two resolutions: a low resolution with cell spacing 
A7?'/2 = 0.2AU'/^ and 1000 grid cells, and a high resolution 
with with A/?i/2 = O.OSAU'/^ and 4000 grid cells. In the ab- 
sence of a planet the low resolution is used, but we switch to 
high resolution when a planet is present. Switching to high 
resolution is achieved by simple linear interpolation, which 
is sufficient to conserve mass and angular momentum at very 
high accuracy. 

The time-step is typically limited by the radial velocity of 
the gas very close to the location of the planet, and evolv- 
ing the system on such a short time-step is unnecessary after 
the planet has opened a gap in the disk. Consequently we 
impose a maximum torque (and therefore a maximum gas ve- 
locity in the radial direction) of |A| < QARHVL^. In addition, 
we make no attempt to model the interaction of the planet 
with the inner edge of the gas disk (which is truncated by the 
stellar magnetosphere), and simply remove planets when they 
migrate to radii a < 0.15AU. The cumulative numerical er- 
rors in the conservation of mass and angular momentum are 
typically < 0.1% over the lifetime of the model. 



2.3. Planetary accretion 

The description of planetary migration above does not al- 
low any material to flow across the gap in the disk induced by 
the presence of a planet. However, numerical simulations of 
the disk-planet interaction show that tidal streams of gas do 
flow across the gap, allowing accretion on to both the planet 
and t he inner disk to persist after the gap has been opened 
(e.g.. lArtvmowicz & Lubowlll996h . The rate of accretion on 
to the planet can be parametrized as the fraction e of the disk 
accretion rate that would be measured in a steady disk in the 
absence of a planet. Simulations show that the efficiency e 
varies strongly with the mass of the planet, and can approach 
(and even exceed) unity f or planets which are marginally able 
to open a gap in the disk (iLubow et al.lll999t ID' Angelo et al.l 
12002). 

FiglT]shows values of e obtained from two-dimensional nu- 
merical simulations of planets embedded in gaseous disks. 
Simulations of low-mass and high-mass planets have been 
normalized by demanding consistency at Mp = IMjup, and fur- 
ther normalized so that the maximum value of the efficiency, 
Cmax, is equal to unity. In addition, we also show a simple 
fittin g function which provides a good fit to the numerical le- 
sults dVeras & Armitagell2004h . The fitting function takes the 
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Fig. 1. — Relative efficiency of accretion from the dislc on to tlie planet, 
plotted as a fun ction of the planet ma ss. The crosses show the numer- 
ical results from ID'Ang elo et al.1 j2002l) . and the filled circles the results 
from Lubow et al. ( 1999). The solid line sho ws t he fitting formula from 
Veras & Armitage (2004) described in Equation 110! The lowest planet mass 
we consider in our models is O.SMjup. 

form 
e(Mp) 



= 1.67 
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+ 0.04. (10) 



More recent numerical studies suggest that, in addition to al- 
lowing accretion on to the planet, tidal streams also permit 
accreti on across the gap from the o uter to the inner disk. Fol- 
lowing iLubow & D'Angelol (l2006h . we define the accretion 
rate across the gap to be 



1 + e 



where the accretion rate on to the planet is 



Mo 



(11) 



(12) 



:e(Mp)Mdisk. 

Operationally, we compute the disk accretion rate expected 
in the absence of a planet, Mdisk as the accretion rate at three 
times the radius of the planet: 



MdLsk = 37rz^(3fl)I](3a) 



(13) 



At each timestep we subtract a mass dM = dt{Mp + Minner) 
from the cell(s) immediately outside the gap, and add the ap- 
propriate fractions of this mass to the p lanet and inner gap 
edge. For consistency with the results of iLubow & D'Angelol 
(.2006) we adopt emax = 0.5, but note that our results are not 
very sensitive to the exact value of e^ax- This procedure does 
not explicitly conserve angular momentum, as mass is moved 
from larger to smaller radii without regard for the different 
angular momenta of these orbits. What s hould happen to this 
"excess" angular momentum is not clear (iLubow & D'Angelol 
|2006|) . but tests show that if all of the excess angular momen- 
tum goes to the planet (the extreme case), the migration rate 
is slowed only at the 1-2% level. This is smaller than many 
of the other uncertainties in the migration model, so we are 
satisfied that our planetary accretion procedure is robust. 

2.4. Code Tests 

In order to test the accuracy of our numerical method, 
we performed a comparison with the semi-analytic results of 
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Fig. 2. — Final planet radii as a function of formation time. The points 
show our numerical results: black and red represent models with a = 0.01 and 
a = 0.003 respectively. Circles denote the model where accretion across the 
planet-induced gap is permitted (emax = 0.5); crosses the model where no gas 
flow s across the gap. The lines show the corresponding curves from Armitage 
(20071) . for both partially (dotted) and fully suppressed (dashed) migration 
models. The small discontinuities in the numeiical data for the "flow" models 
at final radii ~ 3AU occur when the planet triggers disk clearing, as discussed 
in the text. When gas does not flow across the gap this triggering is much 
more dramatic, and no planets survive at radii < 1.5AU. 



lArmitagj ([2007). For this test we set the initial disk mass to 
be Md = 10"' ^Mo ~ O.O316M0, and "form" planets of initial 
mass IMjup at a = 5AU. Figure |2] shows the final planet radii 
as a function of formation time, for disks with a = 0.01 and 
a = 0.003 and models with e,nax = 0.5 and e^ax = (i.e., with 
and without accretion flow across the pla net's orbit). Als o 
shown are the semi-analytic predictions of' Armitagg (l2007h . 
with parameters adjusted to match the disk model used here, 
for both "partially suppressed" and "fully suppress ed" models 
(see also Syer & Clarke 1995: llvanov et al.l ll999l)^ In both 
cases the numerical results agree very well with the predic- 
tions of the fully suppressed migration model, so we are sat- 
isfied that our numerical procedure is accurate. The fact that 
our models appear consistent with the fully suppressed migra- 
tion model is not altogether surprising, because by allowing 
gas to flow across the gap we prevent the pile-up of material 
near the outer gap edge that drives systems into the partially 
suppressed regime. 

We also see from Fig|2]that the inclusion of the direct pho- 
toevaporative wind results in a degree of coupling between 
the formation of planets and the onset of disk clearing. If a 
planet suppresses accretion in the inner disk sufficiently, the 
resulting gap in the disk can become optically thin to ion- 
izing photons. As a result the wind switches to the direct 
regime and clears the disk, preventing further planet migra- 
tion; effectively, the planet triggers disk clearing. This be- 
havio ur is distinct from that considered in previous models 
(e.g., lArmitage et al.ll2002l: lArmitagell2007h . which assumed 
that disk clearing was independent of planet formation and 
migration. When accretion across the gap is significant the 

^ Suppression of Type II migration occurs when the planet mass Mp ex- 
ceeds the local disk mass 47rSc!^, as the planet's inertia then causes it to 
migrate on a time-scale longer than the local viscous time-scale. The degree 
of suppression is characterised by the parameter B = 4nSa^/Mp: for partial 
suppression the migration rate is reduced by a factor B'/^, while in the fully 
suppressed case it is reduced by a factor B. 



consequences of this are not dramatic (as seen in Fig|2|i: the 
location of the small discontinuity in allowed final planet 
radii differs for different disk models, and when we assume 
a spread in disk initial conditions the discontinuity is not re- 
flected in the resulting distribution of planet radii. However, 
when accretion across the gap is suppressed (i.e., for very low 
values of e) we see strong coupling between planet formation 
and disk clearing, with significant implications for the final 
distribution of planet semi-major axes. In this case planets 
can only survive if the outer edge of the gap they induce is 
outside the critical radius where photoevaporation opens a gap 
in the disk (approximately 2AU); at smaller radii, the planet 
is always swept on to the star when the inner disk is cleared. 
Consequently, planets can only survive at radii < 1.5AU if 
there is an accretion flow across the gap. Such planets are 
commonly observed, so we conclude that some accretion flow 
across planetary orbits must occur in real systems. 

2.5. Reference model 

In order to study the effects of migration and accretion on 
the disk and planet properties, we run sets of models in which 
we randomly sample various model parameters. We first de- 
scribe our reference model set, and then a number of further 
model sets in which we explore the effects of varying different 
parameters in the model. 

In the reference model set the parameters of the disk model 
are fixed to the values described above, but we allow for a 
spread in initial disk masses. We use a log-nor mal distribution 
of initial disk masses (see, e.g., [Kennedy & K envon 2009), 
with a mean of log[o((Md}/Mo) = -1.5 and a 3-(t spread of 
0.5dex. In the absence of a planet the mean disk model has an 
initial accretion rate of ~ 4 x 10~^MQyr"', starts to undergo 
photoevaporative clearing after a "lifetime" of ~ 4.2Myr, and 
is completely cleared after ~ 4.7Myr 

Not all stars are observed to possess giant planets at AU 
radii, so we form planets only in a subset of our models: the 
probability of a giant planet forming in each individual model 
is 10% (this value chosen arbitrarily in order to reproduce the 
observed frequency of giant planets). In each model which 
forms a planet, we then allow a single planet of mass Mp to 
"form" at time fp and radius flp. In the reference model set we 
use a constant planet formation radius of Op = 5AU. We as- 
sign the initial planet masses by randomly sampling a distri- 
bution p(Mp) oc 1 /Mp (e.g., Marcy et al. 2008), over the range 
0.5Mjiip < Mp < 5.0Mji,p. Here the lower hmit of 0.5Mjup is 
approximately the minimum gap-opening mass for our disk 
model, and the upper limit is chosen, somewhat arbitrarily, 
near the limit set by the most massive known extra-solar plan- 
ets. 

The time of formation is assigned randomly in the range to 
0.25Myr < fp < fc. The lower limit represents the earliest time 
at which our disk model is physical (see ^2.1l above). and the 
upper limit is the time at which the wind begins to clear the 
gas disk. Following Clarke et al., (.2001.) and .Ruden (.2004i) we 
define this time as 



_ 1 f 3Mi 



2/3 



(14) 



In addition, we limit the maximum value of tp such that the 
planet mass cannot exceed the instantaneous disk mass at tp. 
In practice, however, the wind begins to clear the gas disk 
when the disk mass reaches ~ 5-lOMjup, so this constraint 
only applies to the most massive planets. We ran = 1000 
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M„ > 1.65M, „, 0.1AU < a < 2.5AU 

p Jup' 




a / AU 

Fig. 3. — Cumulative d istributions of p l anet s emi-major axes. The black 
line shows the data from Fisch er & Valentil <2005l) . cut to be complete in both 
mass (Msinf) and semi-major axis. The red line is the coiTesponding distri- 
bution from our reference model. The dotted red line shows the distribution 
from a second random realisation of our reference model, and illustrates the 
typical Poisson errors associated with samples of this size. 

randomly realized models: a total of 93 models formed plan- 
ets. 54 of these planets survived; the remainder were accreted 
on to the central star. 

2.6. Results 

2.6.1. Planet properties 

A first test of this approach is to reproduce the ob- 
served semi-major axis dist ribution o f extra -solar planets. 
We follow the approach of lArmitagd (|2007|) in creating a 
uniformly-selected, complete sample of extra-solar planets, 
us ing the data from the L ick radial velocity survey given 
in iFischer & Valentil (l2005 b"^. These data are complete for 
Doppler velocities K > 30m s"' and orbital periods P < 4yr. 
Consequently, we apply cuts to the sample so that Msin/ > 
1.65Mjup and a < 2.5AU. Our model makes no attempt to 
model the survival of "hot Jupiters" at small radii, so we 
also require that a > 0.1 AU. This results in a final sample 
of 23 extra-solar planets (from 850 host stars). The radial dis- 
tribution of these planets is shown in Figl3] along with the 
corresponding distribution from our reference model (for the 
21 surviving planets with final masses > 1.65Mjup and semi- 
major axes < 2.5AU). The two distributions are qualitatively 
similar, and a Kolmogorov-Smirnov (KS) test fails to reject 
the (null) hypothesis that that the two data sets are drawn 
from the same underlying distribution. (The KS probability 
is 10%, with a KS statistic of D = 0.35.) The apparent lack 
of planets around 1 .5-2AU is a chance fluctuation, and is not 
statistically sign ificant (see also Fig|2]and the associated dis- 
cussion in ^2.4l i. Also shown in Figl3]is the distribution re- 
sulting from a second randomly realized set of 1000 models: 

Where possible, we have updated the data fromlFisc her & V alentil j200S) 
using the improved orbital parameters given in IButleretalJ<2006l) . 



this time 105 planets formed, and 56 survived. The null hy- 
pothesis is again not strongly rejected by a KS test, but the 
KS probability of 2% is somewhat lower than before: this il- 
lustrates the typical Poisson errors associated with data sets 
of this size. A less conservative cut of the data (K > 20m s"', 
P < 5yr, Msin / > 1.2Mjup, 0.1 <a< 3AU) results in slightly 
improved number statistics (33 exoplanets), but quantitatively 
similar distributions of planet semi-major axes. 

As a further test, we confirm that the frequency of plan- 
ets in our models is als o consistent wit h the observed ex- 
oplanet statistics. Gumming et al.l (l2008i) report frequencies 
of 2.1 ±0.7% for planets with Msin; > 2Mjup and O.IAU 

< a < 2AU; 0.6 ± 0.4% for Msin/ > 2Mjup and 2AU <a< 
3AU; 1.3 ±0.5% for IMjup < Msin/ < 2Mjup and O.IAU 

< fl < 2AU; and 0.6 ±0.4% for IMjup < Msin/ < 2Mjiip and 
2AU < a < 3AU^. The corresponding frequencies for our ref- 
erence model are 1.3%, 0.1%, 2.4% & 0.6% (13, 24, 1 & 6 
planets respectively), which show good agreement with the 
data. However, it should be noted that both the absolute fre- 
quency and the mass function of planets are primarily deter- 
mined by their input values, and are essentially free parame- 
ters in our model. We have chosen these in a simple manner 
in order to be consistent with the observed data, but the ra- 
dial distribution of planets is the more powerful test of migra- 
tion physics. In addition, there is a suggestion that our model 
over-produces planets of l-2Mjup at radii < 2AU. Given the 
small numbers of planets this discrepancy is not statistically 
significant, but it highlights the need for better statistics in this 
field. Larger statistical samples of exoplanets will dramati- 
cally increase our ability to discriminate between models, and 
have the potential to provi de precise const raints on theories of 
planet migration (see also lArmitagell2007l) . 

We are therefore confident that our reference model cor- 
rectly reproduces both the frequency and the radial distribu- 
tion of known extra-solar giant planets, at least at the ~ l-cr 
level. The fact that simple mod els can reproduce this distribu- 
tion has been noted before (e.g.. I Armita gel2007l) . but previous 
such studies have used analytic methods, rather than integrat- 
ing each model explicitly as we do. Our method instead en- 
sures that the disk properties (surface density, accretion rate, 
etc.) of every model are known throughout, and this "brute 
force" method allows us to study the effects of migrating plan- 
ets on the disk population in a manner that is n ot possible 
using more typical population s ynthesis methods (llda & LinI 
l2004allbl: iMordasini et al.ll2009allbh . 

2.6.2. Disk properties 

In addition to studying the properties of the migrating plan- 
ets in our model, we are also able to study the evolution of 
observable disk properties over the lifetime of the model. In 
particular, we are able to follow the accretion rate on to the 
star (i.e., at the inner boundary) and the disk fraction as func- 
tions of time, for our compete set of 1000 disk modes. The 
evolution of the accretion rate is shown in Fig|4l the color 
scale denotes the probability of finding an individual disk at 
any given position in the M-t plane^. The accretion rates de- 
cline from a median value of ~ 5 x lO~^M0yr"' at f ~ lO^yr 

^ The quoted unce rtainties here are purely Poisson errors. In the 
ICumming et alj j2008l ) sample of 475 stars, the (completeness-corrected) 
numbers of objects in these four independent bins are 10, 3, 6 & 3 respec- 
tively. 

* Note that this probability is computed as a fraction of the total lifetime of 
our models, and therefore does not account for the duration of the disk-less, 
WTT phase. 
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Fig. 4. — Evolution of accretion rate as a function of time in the reference 
model (upper panel). The color scale shows the probabihty of any individual 
model being found at any given point on the plot. We regard lO^'^MQyr^'as 
an arbitrary sensitivity limit, denoted by the dotted line: pixels immediately 
below this line include all points with accretion rates less than 10"'"MQyr"' , 
and can be regarded as loci of observational upper limits. The green dashed 
line shows the evolution of the m edian disk mode l in the absence of a planet. 
The data points are taken from IHartmann et alj (1993). T he low er panel 
shows the corresponding results for the SCALE model (see j|2.6.3t : the ef- 
fects of increased dispersion in the disk model are clearly seen. 




t / Myr 

Fig. 5. — Disk fraction as a function of time. The black Hues are from 
our reference model; the red lines from the SCALE model. In both cases the 
solid line is the disk fraction (At /Mot; see and the dotted line is the 
fraction of transitional disks (Wt/Moi; see The data points and error bars 
are the observed disk fractions for a number of nearby sta r-forming regions, 
measured primarily by near-infrared excess, compiled by Mamajek (2009|). 
[These data also have significant systematic uncertainties in the derived ages 
(typically itlMyr), but for clarity we have omitted these en'or bars.] The 
dashed lines shows the disk fractions multiplied by 0.8, and represent the 
expected infrared excess fractions once close binaries are taken into account. 

to 2± lO~'*'M0yr"' at f ~ 4Myr, at which point the photoevap- 
orative wind becomes dominant and the accretion rate drops 
precipitously. There is also significant scatter in the evolu- 
tion: the shortest disk lifetime is 2.3Myr, while the longest is 
lO.VMyr. Much of this scatter is due to the intrinsic dispersion 
in our disk model, but the effects of migrating planets increase 
the scatter significantly compared to a disk-only model. The 
accretion rates in the model are somewhat lower than in the 
IHartmann et alJ (11998! ) data, by a factor of ~ 3. However, this 
is within the systematic uncertainties associated with both the 
measured stellar ages and accretion rates, and given the addi- 
tional uncertainties in our understanding of angular momen- 
tum transport in disks this discrepancy is probably not signif- 
icant. In addition, the maximum disk masses in the model are 
~ O.12M0, with a median val ue of 0.0045Mf7), in good ag ree- 
ment with the observations of Andrews & Williamsl (1200 5). It 
has previously been shown that viscous accretion disk models 
are broadly consistent with the observed evolution of accre- 
tion r ates in CTTs (e.g. . IHartmann et anil998l : lArmitage et alj 
HOOI, and our results re-confirm this conclusion. 

Fig|5] shows the decline in the disk fraction as a function 
of time. We define the disk fraction as the fraction of disks 
of a given age that are "normal" viscous disks. Disks con- 
taining migrating planets, and disks which are undergoing 
photoevaporative clearing, are excluded from this definition 
and are instead classed as "transitional" disks (see ^ be- 
low). This definition should be broadly consistent with pre- 
vious observational studies of the disk fractio n in young star- 
formi ng regions (e.g.. .Haisch et al..2001: Sicilia-Aguilar et al.l 
|2006|) . which identify disks through near- or mid-infrared ex- 
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cess emission (which is suppressed or absent in both transi- 
tional disks and WTTs). We see that the disk fraction declines 
from neai-ly 100% at t < 2Myr to almost zero at t > 8Myr, 
with a median disk lifetime of approximately 4Myr This be- 
haviour agr ees qualitatively with observational studies (e.g., 
iHaisch et al. 2001: Mamaj ek 2009). but results in a disk life- 
times that are somewhat too long (although, given the large 
systematic uncertainties in determining the ages of pre-main- 
sequence stars, it is not clear if this discrepancy is signif- 
icant). However, our model set considers only the evolu- 
tion of disks around single stars, and a significant fraction 
of young stars are binary or multiple systems. The evolu- 
tion o f disks in binary systems is complex (e.g..l Monin et al.l 
l2007h . and beyond the scope of this paper, but we note that 
CTT binaries with small separations (< lOAU) show signif- 
icantly reduced infrared excesses compared to single CTTs 
(due to tidal disruption of the inner disk). Disks in binary 
systems ar e interesting in their own r ight (see, e.g., the dis- 
cussion m [Kennedy & Kenyonll2009h . but for our purposes 
these binaries can be regarded as contaminants in photometric 
studies of disk fractions in young clusters. Some nearby star- 
forming regions, such as Taurus- Auriga, have been the subjec t 
of detailed multiplicity surveys (e.g.. iWhite & Ghezll2001h . 
but more distant and more heavily embedded regions have 
not been similarly studied. In addition, recent surveys us- 
ing high-resolution imaging have discovered significant num- 
bers of close binaries even in regions that had previously been 
well studied ([Ireland & Krausll2008l: iKraus et al.ll2008l) . sug- 
gesting that the fraction of stars that exist as binaries with 
separations of ^ 3-3 0AU is not well constrained. Com- 
parisons to field stars jDuquennov & Mavor 1991) , and pre- 
liminary results from new observational studies tKraus et alj 
l2009i) . suggest that the fraction of binaries with such separa- 
tions is approximately 10-20%. We note also that even in the 
youngest clusters the IR excess fractions rarely exceed 80- 
85%, andjire^inconsistent with 100% even at ages < 0.5Myr 
(e.g.. [Haisch et al.l2001HMama iek 2009). The suggestion that 
some young stars are born without disks poses serious ques- 
tions of our understanding of star formation, but a simpler in- 
terpretation is that the majority of these "disk-less" objects are 
in fact binary or multiple systems. If we assume that 20% of 
young stars are binaries which do not show strong infrared ex- 
cesses, then our model agrees very well with the observed data 
(see Fig|5]l. We do not reproduce the lon g tail of the distribu- 
tion observed in some older clusters (e.g.. lLawson et al]|2004l : 
[Sicilia-Aguilar et al. 2006), but realistic dispersion in some of 
the parameters held constant in our models (a, $, H/R), or 
dispersion in the ages of stars in individual star-forming re- 
gions, could easily produce such a tail. Moreove r, some long- 
lived disks are known to be circumbinary (e.g.. i Furlan et alJ 
l2007h . so "contamination" by binaries probably plays a role 
here too. We are therefore satisfied that our reference model 
successfully reproduces both the observed distribution of gi- 
ant extra-solar planets, and the observed properties of disks 
around young pre-main-sequence stars. 

2.6.3. Effects of model parameters 

In order to study the effect of various parameters on the re- 
sults of our modeling, we also ran a number of model sets 
with parameters different from those of the reference model. 
In the reference model all planets form at a fixed radius 
ap = 5AU. In reality we expect giant planets to form over a 
range of radii, so we also considered sets of models where ap 
was assigned randomly in the ranges 2AU< Op < lOAU and 
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TABLE 1 

List of model sets ran. For each model set 1000 individual models were run, 
with planets forming in 10% of the disks. The values listed for the initial 
disk mass are the means of the log-normal distributions sampled, as 
described in the text. Where the other listed parameters were not fixed, they 
were randomly drawn from the distributions described in the text 
[p(Rs) = constant, p(ap) = constant, p(Mp) oc 1 /Mp]. 



5AU< flp < lOAU. In addition, we consider three other vari- 
ant model sets: one with logjo((Md)/MQ) = -1.0: one with 
a = 0.003: and one with the scale radius assigned randomly 
in the range 5AU< Rs < 20AU. (This last variant has the ef- 
fect of creating a dispersion in the characteristic viscous scal- 
ing time fjy.) Lastly, we ran a set of models where planets 
formed with a fixed mass of O.SMjup, rather than the distribu- 
tion of masses used in the reference model. The parameters 
of our model sets are listed in full in Table [T] 

The resulting distributions of planet radii unfortunately do 
not discriminate significantly between the various models. KS 
tests fail to reject any of the models which allow for a range 
of initial planet masses (probabilities in the range 0.5-15%), 
and on this basis none of these models are strongly preferred 
(or disfavoured) with respect to the reference model. The fre- 
quencies of surviving planets vary only weakly between the 
models, and given that this depends primarily on the input 
planet formation frequency this also fails to constrain on the 
model parameters significantly. The low-viscosity (ALPHA) 
and high disk mass (DISKMASS) models are somewhat dis- 
favoured due to their long disk lifetimes (mean disk lifetimes 
of 7.4 & 9.8Myr respectively), but given the intrinsic uncer- 
tainties in our understanding of angular momentum transport 
in disks this is also not especially significant. Similarly, the 
SCALE model, in which the disk scale radius R, varies, is 
weakly preferred over the reference model, due to the in- 
creased dispersion in the resulting disk lifetimes and accre- 
tion rates (see FigHjl. Allowing for ranges in planet formation 
radii does not significantly alter the distribution of planets in 
the "migration zone" (< 3AU), but does give rise to signif- 
icant differences in the distribution of planets at larger radii 
(~ 5-lOAU). When the radial velocity surveys for planets be- 
come complete to larger radii (longer orbital periods), they 
will provide stronger constraints on the radii at which giant 
planets form (see also Armitage 2007). 

The one model which is strongly rejected is that in which 
planets form with a fixed mass of 0.5Mjup (model FIXED- 
MASS). In this model dispersion in planet masses is solely 
due to differences in the accretion history of planets during 
the migration phase, and consequently the final planet masses 
correlate strongly with radius. Planets which spend a longer 
time migrating accrete more gas, so the most massive plan- 
ets are always found at small radii: such a correlation is not 
observed in exoplanet surveys. Moreover, planetary accretion 
is not fast enough to account for the observed range in planet 
masses: the most massive surviving planets have masses of 
~ 1.5Mjup. We are therefore able to reject this model at high 



8 



Alexander & Armitage 



confidence, and conclude that giant planets must enter the 
Type 11 migration regime with a large range of masses. 

2.7. Discussion 

A critical feature of our migration model is the maximum 
time at which planets are allowed to form, tc- Planets which 
form at late times are only able to migrate a limited distance 
before the disk is cleared (see Fig|2|l, so the latest time at 
which planets can form has important implications for the 
resulting distribution of planets, especially close to the for- 
mation radius Op. Variations in tc do not strongly affect the 
distribution of planets in the reference model at radii < 3AU, 
but changing by as little as 10% can result in 3-cr changes 
in the distribution at larger radii. We therefore attach lim- 
ited significance to our predicted distributions of planets in 
the range ^3-lOAU. However, we note that even the forma- 
tion of the observed population of giant planets (in the mi- 
gration zone) requires planets to form at relatively late times, 
>l-2Myr (see Fig|2]i. This result is predicated on the assump- 
tion that protoplanetary disks accrete in a manner consistent 
with our viscous accretion model. Adopting lower values 
of a reduces the efficiency of Type 11 migration somewhat, 
but one cannot adopt arbitrarily low values and still repro- 
duce the observed accretion rates. As long as disks are ac- 
creting viscously at the observed rates, it is difficult to avoid 
the conclusion that giant planets must be able to form late. At 
these times the disk masses in our models are low, < O.OIM0, 
with surface densities at 5-lOAU of < lOg cm~^. This is at 
least an order of magnitude smaller than the s urface density in 
the ca nonical Minimum Mass Solar Nebula ( WeidenschillingI 
Il977 t Hayashi 1981), which is treated as a fiducial value 
in many calculations of planet formation. Forming planets 
in disks with such low surface densities may be challenging 
for m odern theories of planet formation (e.g.. lJohansen et aH 
lIOOTh . 

A further limitation of our models is that we consider only 
the formation of one planet per disk, while in reality many 
planetary systems contain multiple planets. In such systems 
planet-planet interactions can be important, and can modify 
the extra-solar planet distribution after the dispersal of the gas 
disk. Indeed, it seems likely that planet-planet scattering is 
responsi ble for the observed dist r ibution of exoplanet ec cen- 
tricities ("J uric & Tremain3l2008l: IChatteriee et al.ll2008 h. In 
this scenario many planets in the migration zone undergo ad- 
ditional migration after the gas disk is cleared, and this pro- 
cess causes the innermost planet to migrate by an amount that 
depends, on average, on the number of giant planets present 
at the end of the disk lifetime. If this number is typically 
small (2-3) scattering will result in a modest re-mapping of 
our planet distributions to smaller radii, but will not remove 
the need for the distribution of planets at small semi-major 
axes to be primarily established through disk migration. The 
shape of the distribution would also change if the scattering 
properties are not scale-free in radius. Considering these ef- 
fects would be necessary in precision tests of data against the- 
ory, but is not warranted with the limited suitable data samples 
available at present. 

3. TRANSITIONAL DISKS 

Since their discovery by IStrom et"aLl (1 19891) . the so-called 
"transitional disks" have been thought to represent a cru- 
cial step in the evolution of planet-forming disks around 
young stars. These disks are characterized by reduced emis- 
sion in the infrared but longer wavelength emission con- 



sistent with normal CTT disks, and this is generally at- 
tributed to some degree of inner disk clearing. The frac- 
tion of disks which appear tra nsitional is small, typically 
5-10% (Skrutskie et al. 1990: iKenvon & HartmannI 119951: 
[Andrews & Williams 2005), which le ads to the conclusion 
that the transitional ph ase is short-lived (ISimon & Prato 1995t 
IWolk & Walt"erlll996l) . Unfortunately this also means that the 
sample of well-studied transitional disks is small, and this 
small sample shows considerable diversity in disk properties. 
A detailed understanding of this important phase of disk evo- 
lution^ has thus so far remained elusive (see the discussion in 
|Ale2Lander2008). 

Theoretically, a number of different physical processes 
are expected to give rise to disks w ith "transi t ional" 
spectral energy distributions: planets (Rice et al. 20031 
Ouillen et al. 2004), dust evolution (Dulle mond & Domin iM 
2005: Krauss et al. 2007 |), disk clea ring ([Alexander et al.l 
2006b; Chiang & Murrav-Cla vl |2007|) and the presence of 
companions (Jensen & Mathie u|[l997l: [Ireland & Krausll2008l) 
probably all play a role. However, recent observational stud- 
ies, especially those with the Spitzer Space Telescope, have 
led to a dramatic increase in the number of known transitional 
disks, resu lting in the first demographic studies of their prop- 
erties. N aiita et all (l2007b identified a sample of 12 transi- 
tional disks in the Taurus-Auriga cloud, using 5-30/im data 
from the Spitzer spectroscopic survey of Furlan et al. ( 2006). 
All of the 8 single stars in their sample are actively accret- 
ing, but the accretion rates for the transitional disks were 
found to be, on average, an ord er of magnitude low er than for 
CTTs with similar disk masses. Naiita et al.l (l2007h concluded 
that partial inner disk clearing by embedded planets was the 
most probable explanation for this result, but noted that other 
explanations (notably dust settling or gr owth) could not be 
ruled out. By contrast. ICieza et alT(l2008l) identified a sample 
of 26 transitional disks in a number of nearby star-forming 
clouds, using different selection criteria based on observations 
across a wider range in wavelengths. They found that inner 
disk clearing was associated with significant depletion of the 
(outer) disk mass, and concluded that disk evolution (presum- 
ably due to viscosity and/or photoevaporation) was the most 
probable explanation. These contrasting results suggest that 
selection biases still dominate these relatively small samples, 
but also suggest that more than one physical mechanism is 
responsible for the systems that are broadly classed as "tran- 
sitional". This view is further supported by the recent results 
of Salyk et al. (2009), who used observations of CO emission 
lines to divide a sample of 14 transition disks into "cleared" 
and "partially depleted" inner disks. 

A key issue in such studies is how transitional d isks are de- 
fined. A broad definition, such as that used by Najita et al.l 
(l2007h . encompasses settled dust disks, disks where signifi- 
cant grain growth has occurred, disks with inner holes, disks 
with embedded planets, and som e binaries. A mor e strict def- 
inition, such as that proposed by [Alexander! (l2008h . limits the 
sample to objects with partially or fully cleared inner holes, 
and thus only selects objects whose gas disks have undergone 
significant evolution or perturbation. In our models we can 
identify, and distinguish, two different types of "transitional" 
disk: disks with holes or gaps due to embedded planets, and 
disks which are being cleared (through the combined action of 
viscosity, photoe vaporation, and possibly planetary torques). 
iRice et alT (12006) demonstrated that a planet which opens a 
gap in a disk reduces the dust-to-gas ratio in the inner disk, 
resulting in a corresponding suppression of infrared emission. 
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As we only consider planets which are massive enough to 
open such a gap, all planet-bearing disks in our models are 
likely to be identified as transitional. In addition, any disk 
which has evolved to be optically thin in the inner disk at in- 
frared wavelengths but remains optically thick at larger radii 
(i.e., any disk undergoing inside-out clearing) is likely to be 
classed as transitional. Consequently, in our models we iden- 
tify disks as transitional if 

I](1AU) < 10-2g cm-2 AND j:(Rh) > IQ-^g cm'^ 
OR 

the disk contains a planet. 

Note, however, that disks identified in this manner represent 
only a subset of the observed samples of transitional disks, 
which a re g enerally more broadly defined (as discussed above 
and in 33.3b . 

In our models the disks are able to spread to arbitrarily large 
radii, and many of our disks expand to radii of order lOOOAU. 
However, such large disks are not commonly observed (e.g., 
[Andrews & Williams 2007; Andrews et al. 2009), and in re- 
ality the outer edges of disks are likely to be truncated by a 
variety of different physical processes (such as photoevapora- 
tion by non-ionizing radiation, or tidal interactions with other 
stars). This simplification does not have a strong effect on the 
global evolution of our models, but does cause problems in 
defining when disks have been cleared. We define all disks 
which have inner holes larger than some critical radius /?/, to 
be "cleared". Initially we adopt R/, = lOOAU, but we also in- 
vestigate the effect of varying this parameter on our results. 

For clarity we define the following quantities in each set of 
models, which vary as functions of time: 

Mot The total number of stars (always 1000 in our models). 

A^c The number of stars with normal, non-planet-bearing 
disks, analogous to the number of CTTs (or Class II 
sources). 

A't The number of stars with transitional disks (according 
to the above criteria). 

A^P The number of stars with planet-bearing disks. These 
represent a subset of the more broadly-defined transi- 
tional disks. (Note that stars with planets but no disk 
are classified as disk-less.) 

A^v The number of stars with disks cleared to beyond 
Rh, analogous to the number of WTTs (or Class III 
sources). 

By construction, Nc+Nj +Nw = Mot- We note in passing that 
the term transitional is rather misleading in this context, as 
many of our disks evolve from a "transitional", planet-bearing 
phase back into a normal CTT phase before they are finally 
cleared. 

3.1. Number statistics 

By adopting the definitions above, we are able to iden- 
tify a subset of transitional disks in our model sets that is 
broadly consistent with observationally-defined samples of 
transitional disks. The fraction of transitional disks (Mx/Mot) 
in the reference model varies from 0-13%, as seen in Fig|5] 
Varying the value of /?/, has only a small effect on this re- 
sult; for Rh = 30AU the peak transition disk fraction is 10%, 
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Fig. 6. — The fraction of transitional disks in the reference model set 
which contain planets (N^/Nt), plotted as a function of time. From top to 
bottom the three curves are calculated for Ri, = 30. 100 & 300AU respectively: 
selecting transitional disks with larger hole sizes increases the fraction of 
clearing disks. At early times all of the transitional disks contain embedded 
planets. 

while for R/, = 300AU it rises to 17%. The transition fraction 
is initially small, rises to a peak at approximately the median 
disk lifetime, and then declines. The increase in the transi- 
tion disk fraction with time can be understood in terms of the 
increasing migration time-scale (due to increased suppression 
of migration with declining disk mass), and also the increas- 
ing incidence of clearing discs at later times; the decline at late 
times is imposed by overall decrease number in the number of 
disk-bearing systems. However, some aspects of the shape of 
this curve are artefacts of our model: in particular, the lack 
of transition disks at early times is in part due to the fact that 
we do not form any planets at f < 0.25Myr. In addition, we 
assume zero dispersion in the age of our populations, while 
real clusters may have significant age spreads. We therefore 
urge caution when comparing our results to observations of 
very young clusters. 

Fig|6]shows the fraction of transitional disks which possess 
planets {Np/Nj), as a function of time, for three different val- 
ues of Rf, (30, 100 and 300AU). Larger values of R/, naturally 
lower the ratio Np /Nj (because all disks with such large holes 
must be clearing), but the same general trend is seen in all 
the models. At early times, the disk accretion rates are too 
high for photoevaporation to be important, so essentially all 
of the transitional disks with ages < 2Myr are planet-bearing. 
As the population evolves photoevaporation becomes impor- 
tant for a progressively larger fraction of disks, so the num- 
ber of clearing disks increases and the ratio Np/Nj declines. 
Here the time-scale of ~ 2Myr is simply the point at which 
the first of our disks (those with the lowest initial masses) 
reach the low, ^ lO~"^M0yr"' accretion rates where photo- 
evaporation becomes important. The absolute time-scales in 
our model are set by the disk viscosity and initial conditions, 
which are essentially free parameters chosen to match the ob- 
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Fig. 7. — The fraction of disks that are transitional [Nj/iNc+NT); black 
lines], plotted as a function of time (for R;, = lOOAU). The sohd line shows 
the results from the REFERENCE model; the dotted line the SCALE model. 
At late times, when few disks remain, this number approaches unity. Also 
shown for comparison are the "true" transition disk fractions Nj /Wiot , for the 
REFERENCE (solid red line) and SCALE (dashed red hne) models. 

s erve d constraints on disk evolution (as discussed in Section 
^2.5\ . We choose these parameters such that the median disk 
lifetime is ~4Myr, and the time at which the first clearing 
disks appear is determined by the dispersion in disk lifetimes. 
In the reference model set this is set by the dispersion in initial 
disk masses, and the first clearing disks appear after ^ 2Myr; 
in the SCALE model set there is an additional dispersion in 
the viscous scaling time, and the first clearing disks appear 
slightly earlier ('^ l.SMyr). However, unless the dispersion in 
disk lifetimes is very large (of order the median lifetime), it 
is very unlikely that disk clearing will be significant at ages 
< l-2Myr^. These results suggest that younger transitional 
disks may be the most promising candidates for hosting em- 
bedded planets (at least in a statistical sense). They also sug- 
gest that longer wavelength observations, which are sensitive 
to larger hole sizes, should preferentially detect clearing tran- 
sitional disks, and that new facilities (such as the Herschel 
Space Observatory) may discover large numbers of these ob- 
jects. 

FigH plots the fraction of disks that are transitional [i.e., 
A't/(A^c+A't)] as a function of time, for the REFERENCE and 
SCALE models. At early times this "transition fraction" is 
small, but as the total disk fraction declines A^r / (A^c + Ni) 
increases to around 30%, and approaches unity at very late 
times (when only a handful of disks remain, most of which 
are clearing). In addition, the small denominator causes this 
ratio to become highly stochastic at last times. We stress that 
this behavior is not inconsistent with the rapid disk clearin g 
seen in our models (as claimed by, e.g.. ICurrie et al.l [20091) . 
but is instead a natural consequence of a rapid transition in a 

' If, h owever, disk clearing was to begin at high er accretion rates (as sug- 
gested bv lErcolano et alJ2009hHOwen et alJ20O9ll . we may see qualitatively 
different behavior. 
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Fig. 8. — As Fig|4] but for transitional disks only. The green dashed line 
again shows the evolution of the median disk model. The accretion rates for 
accreting transitional disks are suppressed by around an order of magnitude 
relative to the disk population as a whole. 

population of stars with a plausible dispersion in disk lifetime. 
In addition, we note that our models somewhat under-estimate 
the fraction of transitional disks, as we do not consider pro- 
cesses, such as dust settling or terrestrial planet formation, 
that may lead to the ap pearance of "hom ologously depleted" 
transitional disks (e.g.. IWood et al.|[2002i) . Consequently at 
late times, when only a few disks remain, many, if not most 
of them will be transitional. We therefore urge against using 
the ratio Nj / {Nc+ Nj) (or, similarly, Nf/Nc) as a measure of 
the transition disk fraction, especially in clusters older than a 
few Myr (where the disk fraction A^c /Mot is small). Although 
it is observationally less convenient (as it requires an accurate 
census of the number of disk-less stars), the ratio Nj/Ntot is a 
more robust statistical measure of the duration of the transi- 
tion disk phase. 

3.2. Transition disk properties 

Rec ently, [Alexander & Armitag3 ( l2007h and iNaiita et al] 
(l2007h independently proposed that statistical studies of the 
masses and accretion rates of transitional disks could provide 
insight into the physical nature of these systems. Embed- 
ded giant planets are expected to suppress accretion without 
significantly altering the disk mass, while photoevaporative 
clearing requires significant evolution of the entire disk be- 
fore the accretion rate falls to a low enough value for a gap 
to open. However, the theoretical arguments used in these 
papers were highly idealized, and did not consider the time- 
dependence of the planet-disk interaction. Our models allow 
us to study these processes more fully, and make detailed pre- 
dictions about the properties of (some) observable transitional 
disks. 

FiglH] shows the evolution of accretion rates in the refer- 
ence model for disks that are ide ntifie d as transitional. At 
ages < 2Myr (see discussion in ^3.11 above) essentially all 
of the transitional disks are accreting, but at later times there 
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are two coeval populations of transitional disks: accreting, 
planet-bearing disks and non-accreting, clearing disks (a few 
of which also contain planets). In the subset of accreting tran- 
sitional disks the accretion rates show considerable scatter, 
and the median is suppressed by approximately a factor of 10 
with resp ect to the CTT popul ation. This is consistent with the 
results of iNaiita et al.l (|2007|) . and supports their suggestion 
that the transitional disks in their sample possess embedded 
planetSj 

Both lAlexander & Armitage' ('2007") and'N aiita etall ( l2007h 
suggested that the distribution of transitional disks in the 
Mdisk-M plane can be a valuable diagnostic of the properties 
of transitional disks. Fig|9]shows the distribution of accretion 
rates as a function of instantaneous disk mass for the transi- 
tional disks in the REFERENCE and SCALE mode l sets. Also 
plotted are data fr om the studies of iNajita et al.l (l2007h and 
ICieza et al.l (l2008h . As our model has M* = IMq, we only 
plot objects of spectral type Ml and earlier (corresponding 
to stars with > O.5M0). We also omit known binaries 
from the figure, as well as disks whose classification is un- 
certain (those classified "C/T" by Naiita et al. 2007). All but 
one of the remaining stars in the Cieza et al. (2008) sample 
are WTTs with no measured accretion rates; we assign up- 
per limits of 10~'*'MQyr"' to these obiec t s^. Vf e note also 
that the binary statistics of the Cie za et al.l (l2008h sample are 
not well known. We find reasonable agreement between the 
predictions of our models and the data, but find that the REF- 
ERENCE model fails to reproduce the transition disks with the 
highest disk masses. This discrepancy is within the systematic 
errors associated with the observations, but is also in part an 
artefact of the relatively small dispersion in our REFERENCE 
disk model set. The SCALE model set, in which the disk scale 
radius Rs was allowed to vary, shows a larger dispersion in 
the Mdisk-A^ diagram, and provides a better fit to the observed 
transitional disk population. 

All of the accreting (> lO~"M0yr"') transitional disks in 
our models contain planets, while the overwhelming majority 
of the non-accreting transitional disks are undergoing photo- 
evaporative clearing. This is consistent with the predictions 
of pre vious studies (Alexander & Armitage 2007; Naiit a et al.l 
l2007h . but we find that these two populations are not as well- 
separated in the Mjisk-^ plane as previous work has sug- 
gested. In particular, ^ 35% of the planet-bearing transitional 
disks in our models have accretion rates in the range 10"'^- 
10~'°MQyr"'(i.e., are accreting at a rate below the sensitiv- 
ity limit of current observations). This suggests that more 
sensitive observations, with limits of < lO~"M0yr"' will be 
necessary to distinguish cleanly between the populations of 
planet-bearing and clearing transitional disks in this manner 

3.3. Discussion 

These models are the first to attempt to predict the rela- 
tive numbers of different types of transitional disks, but we 
stress that we are not yet capable of modeling all of the 
necessary physical processes simultaneously. In particular, 
our models make no attempt to account for dust settling or 
growth, both of which are known to reduce the disk opacity 
and can give rise to "transitional" spectral energy distributions 
dbulle mond & Dominik 2004, 2005). In essence we consider 
only the subset of transitional disks which have gaps or holes 
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* One object, USco J161420.2-190648, has a measured Ha equivalent 
width of 52A jPreibisch et aU20"02l) but no detected UV-excess, 
it an upper limit of lO^'MQyr"' . 
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Fig. 9. — Distributions of masses and accretion rates for transitional disks. 
The upper panel shows the results from the REFERENCE model set; the lower 
panel the SCALE model set. As in Fig|4] the pixels immediately below the 
dotted line represent upper limits. All of the transition disks with accretion 
rates > 10"' ' Mp)y r~* contain p l anets. Data points are taken from Na iita et al.l 
<2007l circles) and lCieza etal] j20O8l single and double upper limits), with 
the samples cut as described in the text. 

in their gas disks, and consequently our models are not capa- 
ble of spanning the full spectrum of transition disk properties. 
We also note that our "gas-only" definition of "transitional" 
is rather i mprecise; more realistic models of dust-gas cou- 
pling (e.g. . [Alexander &~Armitage,2007) are required to make 
more detailed predictions. As a related point, we urge caution 
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when comparing data to models, particularly with regard to 
selection criteria. All but the most conservative observational 
definitions of "transitional" do not limit themselves to only 
objects with inner holes, and while the objects in these less 
strictly-selected samples are clearly interesting, they cannot 
readily be compared with models of the type presented here. 

We further note that the "planet-bearing" disks identi- 
fied in our models are limited to disks containing gas gi- 
ant planets, which are massive enough to open disk gaps 
and migrate in the Type 11 re gime. Many exo planets of 
lower mass are now known (e.g.,|Udry et al.l| 2007 '), and com- 
parison to the Solar System suggests that terrestrial plan- 
ets can also form on tir ne-scales comparable to observed 
disk lifetimes (see, e.g., iNagasawa et al.l l2007l and refer- 
ences therein). Planets less massive than O.SMjup are un- 
likely to open gaps in their parent gas disks, but the pres- 
ence of even a low-mass planet can cause significant per- 
turba tions to the dust distribution f Paardekooper & Mellemal 
120061) . It therefore seems possible that disks undergoing ter- 
restrial planet formation may also be classed as transitional, 
even in the absence of the more dramatic perturbations to 
the disk structure modeled here. However, the migration of 
low mass planets and planetary cor es is theoretically com- 
plex and not fully understood (e .g., IPapaloizou et al.l 120071: 
iPaardekooper & Papaloizoul l2009l) . and is not constrained by 
current exoplanet data. Consequently any extension of our 
models beyond the Type 11 migration regime would introduce 
significant theoretical uncertainties, and would be of question- 
able benefit in this context. 

In addition, we point out that our models apply only to 
stars of approximately solar mass. Disks around low-mass 
stars and brown d warfs are now commonly observed (e.g. 
IScholzet al.ll2006h . but our knowledge of their evolution is 
limited. Moreover, interpretation of the infrared spectral en- 
ergy distributions of disks around low-mass (M-type) stars is 
fraught with difficulty (Ercolano et al. 2009a), and it is not 
at all clear whether such disks evolve in the same manner as 
their more massive counterparts. Very little is known about 
disk lifetimes and masses in this regime, and the key phys- 
ical processes (angular momentum transport, photoevapora- 
tion, planet formation) are essentially unconstrained. Conse- 
quently, we make no attempt to extrapolate our results to stars 
of lower mass. However, it seems likely that future observa- 
tions will discover a l arge number of transitional disks around 
low-mass stars (e.g.. ISicilia-Aguilar et al.l [20081 : ICurrie et alj 
|2P09), and this represents an interesting avenue for future 
work. 

4. SUMMARY 

We have constructed models of planet migration in evolv- 
ing protoplanetary disks, and used a Monte Carlo approach to 
model the evolution of populations of such disks. The disks 
evolve subject to viscosity and photoevaporation by the cen- 
tral star, and giant planets form and undergo Type II migra- 
tion. Our model successfully reproduces the frequency and 
radial distribution of observed extra-solar planets, and also 
reproduces the observed accretion rates, disk masses and life- 



times of protoplanetary disks. 

The relatively small uniform sample of observed exoplan- 
ets limits the extent to which our models can inform our un- 
derstanding of planet migration, but we are able to draw sev- 
eral interesting conclusions. We find that the addition of the 
"direct" photoevaporative wind results in a degree of cou- 
pling between planet formation and disk clearing, which has 
not been seen in previous models. Consequently some ac- 
cretion flow across the planetary orbit must occur, as other- 
wise is is impossible to "strand" migrating giant planets at 
radii < 1 .5 AU. We also find that planetary accretion during 
the migration phase cannot explain the observed range of ex- 
oplanet masses, and therefore conclude that the planet for- 
mation process must result in a broad range of giant planet 
masses. Lastly, we find that it is only possible for giant plan- 
ets to survive if they form at relatively late times (unless Type 
II migration is dramatically suppressed). At this point in our 
models the disk surface densities at radii of a few AU are low 
(< lOg cm"^), and forming giant planets in such a low-density 
environment may be challenging for current theories of planet 
formation. 

In addition, our models allow us to make a number of pre- 
dictions about the properties and evolution of the so-called 
transitional disks (more precisely, the subset of transitional 
disks with holes or gaps in their gas disks). Our models suc- 
cessfully reproduce the observed transition disk fractions of 
~ 10%, and are also able to explain the accretion rates and 
disk masses of observed samples of transition disks. How- 
ever, we find that the properties of this population evolve sig- 
nificantly with time. We predict the existence of two pop- 
ulations of transitional disks: weakly accreting disks, which 
contain embedded planets, and non-accreting disks which are 
undergoing inside-out clearing. At early times (< 2Myr) es- 
sentially all transitional disks are planet-bearing, but the frac- 
tion of transition disks which possess planets drops with time: 
at late times (> 6Myr), the vast majority of transitional disks 
are being cleared by photoevaporation. Future observations 
will result in much larger samples of transitional disks than 
are currently known, and should allow us to disentangle the 
competing processes of disk evolution and planet formation. 
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